Quantum-based Mechanical Force Realization in Pico-Newton Range 
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We propose mechanical force realization based on flux quantization in the pico-Newton range. 
By controlling the number of flux quantum in a superconducting annulus, a force can be created as 
integer multiples of a constant step. For a 50 nm-thick Nb annulus with the inner and outer radii 
of 5 /^m and 10 /jm, respectively, and the field gradient of 10 T/m the force step is estimated to be 
184 fN. The stability against thermal fluctuations is also addressed. 

PACS numbers: 06.20.fb, 85.25.-j, 84.71.Ba 
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Introduction Due to remarkable improvement in its 
sensitivity, force measurement has become a useful and 
essential probe for leading-edge nano/bio-researchesQ, 
which cover nanoscale imaging, protein folding studies, 
nanoindentation, and many others. The force detection 
limit keeps getting lowered, for example, to an atto- 
Newton (10 -18 N) level in magnetic resonance force mi- 
croscopy. Such an ultimate resolution is enough to read 
a single electron spin. 0, 0, E| 

Unfortunately, however, no direct Si-traceable force re- 
alization has been established even at sub-Newton level, 
which means no standard to compare with measured 
force. Prevailing dead-weight method, which creates 
gravitational force using standard weights, obviously be- 
comes no longer valid below micro- Newton level. Ill At 
the micro-Newton level, an electrostatic force realiza- 
tion has been proposed recently by Pratt et al. Therein, 
electrostatic standard force is created between two coax- 
ial electrodes in maintaining a constant voltage and ex- 
pected to allow a relative uncertainty of parts in 10 4 . Re- 
lated electrical units are traced to their standards based 
on Josephson and quantized Hall effects. At the nano- 
Newton or pico-Newton level, no force realization for 
standard has been suggested despite needs for precision 
measurements of force. It is essential, for instance, to 
testify fundamental forces such as Casimir force, non- 
Newtonian gravitation, etc. 

Also, it casts a striking contrast to the case of electri- 
cal units such as voltage, that (to our best knowledge) 
no attempt has been tried to directly use quantum phe- 
nomena in realizing a mechanical force. Here, we present 
a concept of quantum-based force realization utilizing a 
macroscopic quantum phenomenon known as magnetic 
flux quantization in a superconducting annulus. Mag- 
netic force exerted on flux quanta can be increased or 
decreased by a constant step, which is estimated to be 
sub-pico-Newton level. Moreover, those force-generating 
states are stable and robust against thermal fluctuation 
at liquid helium temperature. 
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FIG. 1: A schematic of force realization based on flux quan- 
tization in superconducting annulus. 



Basic Principles Figure ^ shows a schematic of 
quantum-based force realization. A tens-micron-sized su- 
perconducting annulus is mounted on an ultrasoft micro- 
cantilever. Below a superconducting transition temper- 
ature of the superconducting material, magnetic fluxoid 
through the annulus is quantized, and the resultant mag- 
netic moment has a component with constant steps, the 
number of which depends on the quantum number. The 
step size is determined by fundamental constants such as 
electron charge and the length quantities of the annulus. 

In a calibrated magnetic field gradient, dB/dz, a force 
is created on the superconducting annulus, and therefore 
on the micro-cantilever. Through a designed procedure 
of magnetic field and temperature, we can leave trapped 
flux quanta as many as we want in zero external field, 
while the field gradient is not zero. The displacement 
of the cantilever is monitored by optic interferometer or 
lever and can be used in calibrating its spring constant. 

Magnetic Moment of a Superconducting Annulus For 
the accurate force realization, it is essential first to eval- 
uate the magnetic moment of the superconducting annu- 
lus. 

We consider a superconducting annulus of inner and 
outer radii, a and 6, respectively, and thickness d in a uni- 
form external magnetic field B cxt = B cxt z. In response 
to the external magnetic field or the fluxoid trapped 
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in the hole, currents J are induced in the annulus and 
the field is modified. The net magnetic vector poten- 
tial A = A cxt + A; n( j consists of the external part A ext 
and the induced part Aj n( j. J and A; n( j are related the 
Biot-Savart law 
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where [1q is the magnetic permeability. On the other 
hand, within a superconductor, the current is governed 
by the London equation 
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where 9 is the phase of the superconducting order pa- 
rameter, Al is the London penetration depth, and <I>o = 
h/2\e\ w 2.07 x 10~ 15 weber is the flux quantum. Equa- 
tions and © give an integral equation for J, from 
which one can get the magnetic moment. Given the ge- 
ometry of the superconductor, one can greatly simplify 
the integral equation by following Q . 

Since the current flows only in the thin film (— d/2 < 
z < d/2, d <C Al), it is customary to define a sheet 
current such that 3(x, y, z) w J s h e et(^, y)&{ z ) and to rep- 
resent the physical quantities in the polar coordinates 
(r, (f,z) so that J shcot (a;,y) = J sh oot(r)^, A(x,y,0) = 
[A cxt (r) + A ind (r)]^>, and V9{x,y, 0) = (n/r)(f, where 
the integer n is the number of fluxes trapped in the hole. 
Then the London equation @ reads as 
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where A = \\/d is the effective penetration depth and 
the Biot-Savart law Q is rewritten as 



Ai a d(r) = / dr'Q(r,r') fi J s h cc t(r') , 



where 
Q(r,r') 
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From the London equation ©, one identifies two sources: 
one from the London fluxoid |jj D\{r) = n$ /2irr and 
the other from the external field -D2M = A ex t(r). Ac- 
cordingly, we decompose the induced quantities, J s heet 
and Ajnd, into contributions from D\ and D2, respec- 
tively: J s heet = Ji + J2 and Aj n d = A\ + A<x. Combining 
the two equations J2J and Q, one arrives at the integral 
equations for 3\ and J2 



dr' Q(r,r')^J y {r')+k^J y {r) = -D v (r), (6) 



respectively (v = 1,2) [3 |8{- The nice feature of this 
form is that a < r < b, which facilitates pretty much the 
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FIG. 2: (color online) Total magnetic moment as a function 
of Bext for n trapped fluxoids (solid blue lines), a/6 = 0.5, 
and A/6 = 0.005. The thick (red) line is for the ground state 
configuration. The dashed (red) line is the guide for eyes 
connecting the peaks of ground-state magnetic moment. The 
thick dotted (black) line indicates the experimental procedure 
(see the text). 



numerical solution of the equation. Once J v (y = 1,2) 
are at hand, the magnetic moment m = mi +771,2 is given 
by 



dr r 2 
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We have followed [6( to discretize the integral equa- 
tion © and get J v and m v numerically. 

Figure [21 shows typical behavior of the total magnetic 
moment m of the superconducting annulus as a function 
of Bext ■ It is clearly seen that at a given B ext , the mag- 
netic moment is quantized depending on the number of 
trapped fluxes. 

At a given value of B cxt , the allowed quantized value of 
the magnetic moment corresponds to a metastable state 
of the system (the thin blue lines in Fig-EJ - The magnetic 
moment corresponding to the most stable state (the thick 
red line in Fig. [2Jl is determined by the Gibbs free energy 
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where I = Ii + I2 with I v — J dr J v {r) is the total 
current along the annulus. The ground-state magnetic 
moment is a periodic function of -B cx t with period Bq. 
Figure plots B versus a/b. 

Estimation of Magnetic Moment, Force Step, Force 
Limit To estimate the magnitude of the magnetic mo- 
ment, let us consider an Nb annulus of inner and outer 
radii, a = 5 fim and 6=10 /im, respectively, with thick- 
ness of 50 nm. The Nb material was chosen because 
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FIG. 3: Bo as a function of a/b for different values of A/b. 



its superconducting transition temperature, 9.25 K,[9j is 
readily accessible in cryogenics and its thin film fabri- 
cation is well established to have provided good quality 
of thin- film circuits such as SQUID. Its London pene- 
tration depth Xl and coherence length £ are 50 nm and 
39 nm, respectively. The choice of radii is reasonable 
considering the typical dimension of microfabricated can- 
tilevers with a low spring constant (k = 10~ 4 to 1CP 5 
N/m). ;Hj These parameters fall on the case of a/b = 0.5 
and A/6 = 0.005. 

Then, a magnetic moment step, mQ, by adding a 
single flux quantum, is numerically estimated as 1.116 
TOf,, where the magnetic moment unit mb = <&ob/fio is 
1.65 x 10~ 14 A-m 2 . The force on the superconducting 
annulus is given by 



F 



dB P 



dz 
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and assuming a modest field gradient, dB ext /dz, of 10 
T/m,^J a corresponding force step is estimated to be 
1.84 x 10~ 13 N or 184 fN, which causes ~2nm static 
displacement of a cantilever with k = 10~ 4 N/m. In 
principle, nano- and pico-Newton level can be accessed 
by stacking these steps. 

Practically, a maximum possible force, or the maxi- 
mum number of steps, is limited by the critical current 
of the superconducting material, at which its supercon- 
ducting state cannot be sustained. For the Nb annu- 
lus with above geometry, the force limit is roughly es- 
timated to be ~ 40 pN from the approximate equation 
F m ax — I c K[(a + b)/2] 2 (dB/dz) and the critical current 
density of Nb, j c = 9 x 10 10 A/m 2 at T = 6 K.\V2} 

Force Realization Procedure Realization of only a 
quantum-induced force requires a specially designed pro- 
cedure, because the magnetic moment of a superconduct- 
ing annulus includes a field-induced component linearly 
proportional B ext as seen in Fig. 2. The magnetic field at 
the annulus from the z-gradient magnet, the background, 
and so on is not simple to null satisfactorily, consider- 
ing the small value of Bq (~ 0.13 Oe for the parameters 



above), and would generate a non-zero force offset even 
in zero flux quantum state. Therefore, we need to extract 
selectively the flux quantum contribution. 

In Fig. 2, the difference of magnetic moment, for ex- 
ample, between n — and n — 3 states is always con- 
stant, three-times as large as tuq at any fixed B extl and 
a corresponding change can be realized by transition be- 
tween them in principle. Uncontrolled transitons, such 
as quantum or thermal tunneling (to be discussed later), 
are practically negligible at operation temperatures, ~ 4 
K, due to high energy barriers between metastable states. 

We suggest a force-realization and cantilever- 
calibration procedure as follows: (1) After measuring a 
zero-reference position of a cantilever in n = rii state at 
Bext = and T <C T c , 20] turn on a uniform external 
field up to a target value roughly where n = rif = n^+An 
is a global minimum of Gibbs free energy. But, the 
system still stays at n — rij state due to energy barriers. 
(2) Then, increase temperature above T c and decrease 
it back to T -C T c , which would send the system to the 
global minimum state, n — rif. (3) Decrease and turn off 
the uniform external field and measure the displacement 
of the cantilever. The net change of magnetic moment, 
An x mQ, is only related to flux quanta and so is that 
of a force. From the number change of flux quantum 
(An), annulus dimension, penetration depth, and field 
gradient, the change in force can be evaluated and used 
in determining the spring constant of the cantilever in 
combination with its displacement. 

The starting condition B ext = can be lifted; our pro- 
cedure is still valid even for B ext ^ as long as it is the 
same at the beginning and the end. This is a great ben- 
efit because background magnetic field is hard to shield 
and usually disturbs delicate magnetic experiments. 

In-Situ Calibration of Field Gradient For precise 
force realization, one has to be able to calibrate the mag- 
netic field gradient dB cxt /dz as accurate as possible. Ex- 
perimentally, it is tricky to know the exact value at the 
local position of the annulus on a micro-cantilever, un- 
less the field gradient is strictly uniform in a large volume 
and well pre-calibrated in micro-scale. 

Here we argue that we can calibrate it in-situ by 
means of the characteristic properties of the magnetic 
moment of the superconducting annulus. The magneti- 
zation m(B ext ) in Eq. can be written into the form 



TO 2(-Bcxt) 

m = mi(n) + xB cxt , X= — s > ( 10 ) 

-Dcxt 

where it should be noticed that mi(re) and m2(B cxt ) are 
linearly proportional to n and -B C xt, respectively. The 
field is expanded as 



Bex±\Z) — -Dcxt.O H ; — z ■ 

dz 



(11) 



From Eqs. JSJ, (fTUfl . and (fTTfl . the force F(z) at the po- 
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sition z is given by 

F(z) = [r ni (l)n + xB«t,o] ^ + X f ^ V *• (12) 

az \ dz J 

Then it follows that 

Alj = lAfc = X (dB^y 
uj 2 k 2k, \ dz J ' 1 J 

where u> (fc) is the natural vibration frequency (Hook's 
coefficient) of the cantilever and Aw (Afc) is the shift in w 
(k) due to the field gradient. Let Az = z cq (n) — z cq (n = 
0), i.e., the shift in the equilibrium position z eq for non- 
zero flux n. Since the force F(z eq ) at z = z cq in Eq. 112|) 
should be at equilibrium with (k + Ak)z cq , i.e., 

mi (l)„^E* =kAz, (14) 

one arrives at the relation 

dB ex t 2n mi(l) Aw 

= A ■ ( 15 ) 

az Az x w 

We now note that in Eq. (TSJ, the field gradient dB cxt /dz 
has been expressed only in known quantities of quantities 
that can be determined to high precision: The relative 
frequency shift Aui/uj is typically probed with high pre- 
cision. As shown in Fig.^J the ratio mi(l)/x in Eq. I|15() 
does not depend strongly on the geometry and can be 
easily determined accurately. We therefore conclude that 
Eq. I|15|) allows us to calibrate the field gradient in-situ 
very precisely. 

For k — 1CP 4 N/m and other parameters adopted 
above for a force step estimation, a relative frequency 
shift is estimated to be 2 x 10 -3 from the Eq. I|13(l . and 
becomes smaller for lower spring constant. The relative 
shift is much larger than the sharpness of the resonance 
peak or 1/Q, with a typical quality factor Q = 10 4 ~ 10 5 
for single-crystalline Si cantilevers. |ic| and can be mea- 
sured to high precision. 

Discussion The static displacement Az of the can- 
tilever can be precisely measured using, for example, an 
optic interferometer, which is known to have the rms 
noise less than 0.01 nm in 1 kHz bandwidth at low 
frequencies. 0| However, thermal vibration amplitude, 
Azthi may not be small in comparison with a static dis- 
placement, Azq, due to a single flux quantum. The ratio, 
Az t h/ Azq, is given by y/kksT '/ Fq, where ks and Fq are 
the Boltzmann constant and the force step due to a flux 
quantum, respectively, and is about 0.41 for k = 10~ 4 
N/m and T = 4.2 K, becoming smaller for decreasing 
k and T. But, it may not cause a serious problem be- 
cause the thermal vibration happens most dominantly at 
a resonance frequency (~ a few kHz), especially for high- 
er cantilever, not disturbing the measurement of static 
component. In addition, for large multiples of Fq the 
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FIG. 4: mi(l)/% as a function of a/b for different values of 
A/b. 

relative error due to interferometer and thermal noises 
gets proportionally smaller. 

Instead of a superconducting annulus, a superconduct- 
ing quantum interference device (SQUID) with leads 
could be mounted on the micro-cantilever. Then, instead 
of thermal cycling, a bias current through SQUID can 
be used to permit additional flux quanta, with a merit 
that the whole procedure is done at a fixed temperature. 
Moreover, each entrance of quantum can be monitored in 
real time from the current- voltage characteristics. For its 
benefit to be used, an accurate calculation of magnetic 
moment for the SQUID case is awaited. 

Finally, we briefly remark on the effects of thermal 
fluctuations. In principle, the fluxoids trapped in the an- 
nulus may undergo relaxation processes and their num- 
ber changes. At our working temperature (~ 4 K), the 
relaxation is dominated by the thermal activation, and 
governed by the law 70 exp(— AU/ksT), where AU is the 
energy barrier against the relevant process. For very nar- 
row annulus with b — a <g! A, the underlying mechanism 
is the phase slip 0, 0, 0, llTll. a nd the rate is order of 
several transitions per second |lil Il9| . For a wide annu- 
lus with 6-a> A, the relaxation comes from crossing of 
the vortices across the annulus of superconducting thin 
film. Unlike the narrow limit, there is no detailed study 
of flux-state transitions in this case, neither experimental 
nor theoretical. Naively, the energy barrier can be esti- 
mated by the condensation energy cost inside the vortex 
core AU ~ Hj;£, 2 d/8ir, where H c is the thermodynamic 
critical field. For Nb, H c ~ 0.1 T, and AU ~ 10 5 K.0 
It implies that the effects of the flux relaxation can be 
safely ignored in our scheme using wide enough annulus. 

Conclusion In summary, we suggested pico-Newton 
force realization based on magnetic flux quantum in a 
superconducting annulus. A stepwise force can be gener- 
ated from the annulus on a micro-cantilever in magnetic 
field gradient and the step size depends only on funda- 
mental constants and the length dimensions of the an- 
nulus beside the field gradient. For a 50 nm-thick Nb 
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annulus with inner and outer radii of 5 /Ltm and 10 /im, 
respectively, a force can be created up to ~ 40 pN by a 
step of 184 fN, assuming the field gradient of 10 T/m. 
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